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^ ' Abstract 

' We consider the two-matrix model with potentials whose derivative are arbitrary rational function 

, of fixed pole structure and the support of the spectra of the matrices are union of intervals (hard- 
edges). We derive an explicit formula for the planar limit of the free energy and we derive 

I ' a calculus which allows to compute derivatives of arbitrarily high order by extending classical 

^ . Rauch's variational formulae. The four-points correlation functions are explicitly worked out. 
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1 Introduction 

We consider a matrix model consisting of pairs of Hermitean matrices of size with an (unnormalized) 
probability density of the form 



d/i(Mi,M2) = dMidM2exp 



(1-1) 

ZN{Vi,V2,t) := j dix, t:=Nh. (1-2) 



Here the potentials Vi, V2 are required to have rational derivative; more explicitly we will set 

,a (2^) 

a 

K=i ^ K=i^^ 

V2{y) := V2,oo{y) + YV2M 

a 

V2My) ■■= E ; ^2,«(2/) := E jT^Tj - ^o,« Hy - Pa) (1-4) 

The logarithmic terms in the measure correspond to powers of determinants. Formally the model is well 
defined for arbitrary potentials with complex coefficients, provided that we constrain the spectrum to belong 
to certain contours in the complex plane along the lines explained in [5]. In this case, however, the matrices 
Mi are no longer Hermitean but only normal (i.e. commuting with their Hermitean-adjoint). 
If we insist on a bona fide Hermitean model we should impose that Vi are real functions, bounded from below 
on the real axis. 

In addition to these data we impose that the spectrum contains segments with extrema {Xi} for the first 
matrix and {Yj} for the second matrix (hard— edges of the spectra): in the case of Hermitean matrices then 
we would be restricting the support of the spectra to some arbitrary union of intervals. 
It is used as a working hypotheses that the following limit exists 

HVi,V2,t) := lim ^In^Tv , (1-5) 

TV— »oo iV^ 

where t = Nh is kept fixed in the limit process. 

This model has been analyzed in two papers [12] and [4] from two opposite points of view: in [12] were 
derived the formal properties of the spectral curve and the loop equations in the large N limit, whereas in [4] 
were considered the properties of the associated biorthogonal polynomials and the differentials equations they 
satisfy for finite A^, together with certain Riemann-Hilbert data. 
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The loop-equations show that in the planar limit the resolvents of the two matrices 

W{x)= lim ^/ ^) , W(y)= lim ^/ ^) (1-6) 

^ ' N-.00 N \ x - Ml / ' N^oo N\y-M2/ ^ ' 

satisfy an algebraic equation if we replace y = Y{x) := W{x) — V({x). This means that there is a rational 
expression that defines a (singular) curve in S ^ x -hereby referred to as spectral curve- 

E{x,y)=0. (1-7) 

and that the cuts of the branched covers x : S — > P^ and y : S — s- P^ describe the support of the asymptotic 
density of eigenvalues, and the jumps across these cuts describe the densities themselves. 

From the finite N analysis the spectral curve [5] arises naturally in conjunction with the ODE satisfied by 
the associated biorthogonal polynomials; indeed any S2 consecutive biorthogonal polynomials (were S2 is the 
total degree of the rational function V2(y)) satisfy a (s2 + 1) system of first order ODEs, namely an equation 
of the form 

d,^N{x)=DN{x)^N (1-8) 

and the spectral curve is nothing but Ej\r{x,y) = det{yl — D]\f{x)) = 0. While clearly certain properties 
are valid only for finite or in the infinite limit, certain other properties can be read off both regimes: for 
instance it can be seen [4] that at the hard-edges the matrix Dn{x) has simple poles with nilpotent rank-one 
residue. This implies certain local structure of the spectral curve y{x) above these points. 

In an algebro-geometric approach the functions x,y themselves are meromorphic functions on the spectral 
curve S with specified pole structure and specified singular part near the poles. The loop equations also 
provide a first-order overdetermined set of compatible equations for the free energy; these however are not 
sufficient to uniquely determine the partition function because the polar data of the functions x,y need to 
be supplemented by extra parameters. This is a purely algebro-geometric consideration but they also can be 
heuristically justified along the lines of [7]. It turns out that the extra unspecified parameters can be taken as 
the contour-integrals 

:= ® ydx , (1-9) 

over a maximal set of "independent" non-intersecting contours. The reader with some background in algebraic 
geometry will recognize that there are g =genus(5]) such contours'^. These parameters are often called "filling 
fractions" and in principle they be should uniquely determined by the potentials; the loop equations cannot 
determine the filling fraction but can determine the variations of the Free energy w.r.t. them. This way one 
obtains an extended set of (still compatible) PDEs for T in terms of the full moduli of the algebro-geometric 
problem: we call this function the non-equilibrium free energy. In this situation one can actually integrate 
the PDEs and provide a formula for the planar limit, T. 

Note that the actual free energy of the model is obtained by expressing the filling fraction implicitly as functions 
of the potentials via the equations 

9e,-F(l^i,l^2,t,e) = 0. (1-10) 
Implicit solution yields e = e(Vi, V2, t); the resulting function 

G{VuV2,t) := HyuV2,t,e{yuV2,t)) , (1-11) 



^More appropriately one should consider only the imaginary parts of these integrals over the full homology of the curve. 
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will be called the equilibrium free energy. The distinction is important when computing the higher order 
derivatives of Q, inasmuch as they differ by the higher order derivatives of T by virtue of the chain-rule; indeed 
while 



5Vi{x) 5Vi{x) 



=1(^1,^2,*) 



(1-12) 



(since d^^.J^ = 0) for the second and higher variations the equations do differ, for example 



5^G 



5Vi{x)5Vi{x') 



+ 



5Vi{x)5Vi{x') j^^5Vi{x)5Vi{x) 



=1(^1,^2,*) 



(1-13) 



We will provide simple formulas for both Q,!F. 



The main approach of this paper is similar to that of [1, 2, 13], namely that of ascertain the algebro- 
geometric data in a convenient abstract formulation which provides an explicit formula for T. Once this 
step is accomplished we also build the formalism for the "calculus" that allows to compute arbitrarily high 
order partial derivatives; we recall that the derivatives of represent higher order correlators of the spectral 
invariants of the model in this planar limit and also, the coefficients of their expansion in the parameters of 
the potentials can be related to enumerative problems of polyvalent fat-graph on the sphere. This calculus 
relies on an extension of Rauch's variational formulae to higher order variations (usual Rauch's formulae are 
used to describe first order variations). 

We want to mention an interesting byproduct of the formalism developed for our calculus: indeed, with 
minor modifications mainly in the notation, the same calculus can be applied to computing variation of 
arbitrarily high order of the "tau" function of the universal Whitham hierarchy [19, 22] which in turn is of 
relevance for the Seiberg-Witten model. This application is developed in the Appendix. 

In this application in fact we obtain an extension of the Whitham hierarchy, because on of the two 
primary differentials liX, dY must have some of its simple poles at the points where the other has some of its 
simple zeroes: this requirement follows from the structure of the spectral curve above the hard-edge points. 
In other words in this moduli space the differentials have also some poles at non-marked points on the spectral 
curve. 



1.1 Bergman kernel 

We recall the definition of the Bergman kernel^ a classical object in complex geometry which can be represented 
in terms of prime forms and Theta functions. In fact we will not need any such sophistication because we 
are going to use only its fundamental properties (that uniquely determine it). The Bergman kernel J7(C,C') 
(where C,C' denote here and in the following abstract points on the curve) is a bi-differential on Sg x with 
the properties 

Symmetry: f^(C, CO = ^^(C', C) (1-14) 
Normalization: l 17(C,C') = (1-15) 

*Our use of the term "Bergman kernel" is slightly unconventional, since more commonly the Bergman kernel is a reproducing kernel 
in the space of holomorphic one-forms. The kernel that we here name "Bergman" is sometimes referred to as the "fundamental 
symmetric bi-diffcrential" . We borrow the (ab)use of the name "Bergman" from [15]. 



4 



(f) il(C, CO = 2i7ru;j((;^) = the holomorphic normalized Abelian differential . (1-16) 

It is holomorphic everywhere on x Sg\A, and it has a double pole on the diagonal A := {(, = C'}: namely, 
if z{() is any coordinate, we have 



^ ^ + Isb{C) + o{z{c)-z{C)) 



dz(C)dz(C') , 



(1-17) 



XziO-ziC'W 6 

where the very important quantity Sb{0 is the " Bergman projective connection" (it transforms like the 
Schwartzian derivative under changes of coordinates). 

It follows also from the general theory that any normalized Abelian differential of the third kind with simple 
poles at two points z- and 2;+ with residues respectively ±1 is obtained from the Bergman kernel as 



For later purposes we introduce the dual Bergman kernel defined by 

_ 9 

fi(C, CO := CO - 27r. ^ u;,(C)^fc(CO(B~'),fc 



j,k=i 



where B is the matrix of 6-periods 



(1-18) 



(1-19) 



(1-20) 



In fact J7 is conceptually no different from Q, being just normalized so that <f^ = 0. We keep the distinction 



only for later practical purposes. 



1.1.1 Prime form 

For the sake of completeness we recall here that the definition of the prime form £'(C,CO- 
Definition 1.1 The prime form E{(,(') 's the (-1/2, -1/2) bi-differential on T,g x Sg 



£^(C,CO 



e 



(u(C)-u(CO) 



/ir=i(C)^r.i(CO 



a 









u=0 



(1-21) 
(1-22) 



IS a 



where are the normalized Abelian holomorphic differentials, u is the corresponding Abel map and 
half-integer odd characteristic (the prime form does not depend on which one). 

Then the relation with the Bergman kernel is the following 

9_ 

u^kiCHiC) (1-23) 



f7(C, CO =dcdc' In i?(C, CO = E ^«.^u,lne 

k,j=i 



a 




.P. 





"(C)-u(C') 
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Remark 1.1 In genus zero -of course- there are no theta functions: however there is a Bergman kernel with the same 
properties, given simply by (using the standard coordinate on the complex plane) 




2 Setting and notations 

We extend the setting of the paper [1]^ and we will work with the following data: a (smooth) curve Sg of 
genus g with 2 + K + L distinct marked points oox,^'i, • • • ,Pk, ooy, Qi, - ■ ■ ,qL and two functions X and Y 
with the following pole structure; 

1. The function X has the following divisor of poles 



3. The differential dX vanishes (simply) at the (non-marked) points {^^j and viceversa the differential dY 
vanishes (simply) at the points {rji}. 

All the points entering he above formulae are assumed to be pairwise distinct. The points of the pole divisors 
which are not marked (the rji) will be called "hard-edge" . As hinted at in the introduction these requirement 
follow from either the loop equations [12] or the exact form of the spectral curve [4]: the points Qa := X(qQ,) 
and Xj := X(^j) are the positions of the poles of the derivatives of the potential V-[{X) and the hard-edges 
in the X-plane (and conversely for Y): the fact that the ODE for the biorthogonal polynomials has simple 
poles with nilpotent, rank-one residue at the points Xj, j = 1, . . ., implies that the differential dX vanishes 
at one of the points above Xj, at which the eigenvalue Y has a simple pole. 
Under these assumptions we can write the following asymptotic expansions 




(2-1) 



2. The function Y has the following divisor of poles 




(2-2) 



Y 




K=l 



near ^j, (here Xj 



near pj, (here Q, 



a 




^The functions that there were denoted with P,Q are here denoted with Y,X. 
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X 



X 



/ -25, 



+ near (here Yj := Y{7]j)) 

j=o ^ 

-^2,00+1 ^ 

^^J.ooY'' ^ — + ©(Y ^) near ooy 



J=i 



The asymptotics above imply immediately that there exist two rational functions which we denote by V{ and 
V2 such that 

(2-4) 



Y - VliX) + I) dX , (x- V^iY) + dY 



are holomorphic differentials in the vicinity of the points {co^x., Qa, ct > 1} and {ooy, Pa, ct > 1}, respectively. 
For later reference we spell out these functions 



Vi{x) := Vl^oojx) + (yi,aix) - Uo^aln{x - Qa)) 

a 

K=l K=-\ 



UK. 



- K{x - Q„)^ 



(2-5) 



V2{y) := V2,oo{y) + J2{V2,a{y)-Vo,aHy-Pa)) 



V-, 



J=l 



{ J{y - Pa 



(2-6) 



A local set of coordinates for the moduli space of these data is provided by the coefficients {uK,a,vj^a,t : 
a = 00,1,2,...}, the position of the poles {Qa, Pa}a=i,...: the position of the hard-edge divisors {Xj,Yj} 
together with the so-called filling fractions 



YdX. 



(2-7) 



2.1 Planar limit of the free energy 

The planar limit of the free energy is defined by the following set of compatible equations 



du^ ,T = Uk,o := res X^YdX 

1 1 

duK..^ = UK,a ■■= res — — — -^YdX 
A g„ (X - Qa) 

/■OOx 

duo^^ = Uo,a:=-h YdX 

1 -^"^ 

dx T = Rj := -resY^dX 
2 



a«.^ = -(K,Jx)-^ 



YdX 



dvj.J" = Vjfl := -\ res Y^XdY 

dv,^^ = VK,a ■■= — 7 res— — -yXdY 

J Po (Y - Pay 

/•ooy 

d.o.^=Vo,a:=-h XdY 
1 

dY.T = Si := -resX^dY 

^ 2 r,. 



9p„^ = res y' (Y 



vo,, 



(y-Pa 



XdY 
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^tT = ^^■.= 4- YdX-J2vo,a = 4- XdY-^no,« 

In these formulas the symbol ^f-stands for the regularized integral obtained by subtraction of the singular part 
in the local parameter as follows: 

(i) at oox (ooy) the local parameter is 2; = (z = Y~^); 

(ii) at qa [pa) the local parameter is = X — {z^ = Y — Pa). 

The regularization is then defined as follows: if z is any of the above local parameters then 



:= lim y w - /(e) , 



(2-9) 



where f{z) is defined as the antiderivative (without constant) of the singular part of ^ as a function of z 
(near z = 0). For example 

Example 2.1 The regularized integral according to the definition is 

YdX — lim lim / YdX + \/i.^(X(i?))- I t- Viio.c I ln(X(i?))-"l/i,„(X(e)). (2-10) 



lim lim / YdX + \/i.^(X(i?))- (t-V iio,a I ln(X(i?))-\/i.„(X(e)). 



The two expressions for fi in (2-8) are proven to be equivalent (thus showing the symmetry in the roles of 
X and Y) by integration by parts, paying attention at the definition of the regularization (which involves as 
local parameters X^^ and Y^^ at the two different poles); indeed we have 

YdX = Jm^ ^ YdX - 1^1,00 (X(e)) + {t + J^^^o.^) lnX(e)j = (2-11) 
= Jm^ XdY + X(e)Y(e)-X(p)Y(p) + Fi,oo(X(e)) + (t + ^ uo,a) In X(e) j = (2-12) 

roox 

= _J_ XdY-X(p)Y(p)-(t + ^uo,„) (2-13) 



together with a similar formula for the symmetric expression 



-f YdX = X(p)Y(p) +\t + Y^ vo,(, I - -T XdY . (2-14) 

J ooy V /3 / '^"^ 



Combining the two one has 

/■OOx /"OOy 



M = -/- YdX - ^o,/3 = -f XdY - J2 ^o,a . (2-15) 

/3 



JOOy •'OOx 



In full generality, given any meromorphic differential and local parameters around its poles one can give 
completely explicit formulae for its regularized integrals (see App. D). In our specific setting we give explicit 
formulae of the previous regularized integrals in terms of canonical differentials of the third kind in App. A. 
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We make also the important remark that in order for the above 
formulae to make sense we must perform some surgery on the 
surface by cutting it along a choice of the a, 6-cycles and by per- 
forming some mutually non-intersecting cuts between the poles 
with nonzero residues of the differential YdX. We achieve this 
goal by choosing some segments on the surface joining oox to 
oop, oox to Qa and ooy to p^,. The result of this dissection 
is a simply connected domain where X, Y are meromorphic 
functions and where the regularizations involving logarithms 
are defined by taking the principal determination. 

The compatibility of equations (2-8) for T can be shown 
by taking the cross-derivatives. We now briefly recall, for the 
reader's sake, how to compute them since much of the for- 
malism is needed in the following. The main tool is the previ- 
ously defined Bergman kernel (Sect. 1.1) providing an effective 
way of writing formulas for first, second and third-kind normal- 
ized differentials on the Riemann-surface. This is needed when 
computing the cross derivatives of the free energy since the 
differentials 9YdX and 9XdY (here d is any variation of the 
coordinates) can be identified with certain canonical differen- 
tials. 



bi 




Figure 1 A visualization of an example of the 
dissection mentioned in the text for a genus-2 
curve. 



In order to sketch how, let us first recall the thermodynamic identity 

(9Y)xdX = -(aX)YdY (2-16) 

where the subscript denotes the local coordinate to be kept fixed under variation. As an example of the use 
of (2-16) in identifying the various differentials we consider a derivative 9^^. From the defining relations for 
the coordinates (2-3) we see that 

i'x.M^ - 1 j":^^ ^ "'^"'^ (2-17) 

I C(l)aX near qa 

has a pole of order K at oox without residue. In order what kind of singularity it has at ooy we use (2-16) 
followed by (2-3) 

(8„„Y)xdX = -(8„X)vdY=|«(^_;;>''Y . (2-18) 

I C'(l)dY near pa 

Therefore the differential (9M^^Y)xdX has only a pole at oox and no residues: moreover it follows by 
differentiation of (2-7) that this differential is also normalized (i.e. with vanishing a-cycles), which is sufficient 
to uniquely specify it. It is then an exercise using the properties of ft to see that 

X^ 

(du^ Y)xdX = - res —0 . (2-19) 
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Following the same logic and similar reasoning one can prove the following formulae 



First kind 



(a,Y)xdX = ., = J-£o 



(2-20) 



Second kind 



{du^ ^ Y)xdX = - res —rr^ ='■ ^K,oo 

OOx K 

{du^ „ Y)xdX = res(X - =: u:k,» 

(9x,Y)xdX = resYf]=:a;x, 



(9Q„Y)xdX = res y{,(X 



(X - g„) 



(2-21) 



(^i;jooY)xdX = res —VL =: uj^ 

OOy J ' 



(9,,„Y)xdX = -res(Y - P,)--^J7 = 
(9y^.Y)xdY = - 

(ap„Y)xdX = -res(y2'„(Y) 



res XO =: ujy- 



vo,, 



{Y-Pa 



Third kind 



(9,„,„Y)xdX 
(9,„^Y)xdX 
(9tY)xdX = 



OOx 



OOy 

0, =: 6Jn 



=: U!q 



(2-22) 



The only formulae above that need some further explanations are the ones for the derivatives w.r.t. Xj (or 
similarly Yj); from the asymptotic behavior (2-3) in the local parameter z = y^X — Xj we have 



9x,YdX 



+ 



2zdz 



dz+0{l) 



+ 0(1) 



res Yrj 

6 



(2-23) 

This proves that if 5 = dx^ then the differential has a double pole at without residues: similar reasoning 
at the other singularities and for the a-cycles of the differential force it to be equal to the formula above in 
(2-21). 

Remark 2.1 As explained in the introduction, we are also interested to the restriction of T to the subvariety of the moduli 
space defined by 

d.^T{VuV2,t,e[V^,V2,t)) = Q , j = l,...,g. (2-24) 

Since on this subvariety the differentials XdY, YdX have identically vanishing 6-periods, the formulas for the constrained 
derivatives that substitute (2-21, 2-22)® are the same with O (1-19) replacing Q. 

Before writing the cross derivatives in a way which is symmetric in X, Y, we introduce some useful notation: 
all the differentials (2-20, 2-21, 2-22) are obtained by applying a suitable integral operator to one variable of 

®The equations (2-20) do not make sense on the subvariety since e are not independent coordinates any longer. 
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the Bergman kernel Q according to the following table of translation 



d 



duK,oo 

d 
d 

duK,a 

d 



UK,a '■ = 



1 



roox 

1 



X 



K 



2iKTT L (X - g«)^ 



d_ 

di 
d_ 

dej 



/■OOx 



d 



d 
d 
d 



Vj,oo := 
Vo,a := 
Vj,„ := 



1 



2iJ'K 

Jpc 

1 



OOy 



1 



2ijTT (Y-p^y 
' ^x 



T 



I OOy 



2m 



(2-25) 



All the differentials (2-20, 2-21, 2-22) are obtained by applying the corresponding integral operator in (2-25) 
to the Bergman bidifferential fi. 

In order to write the cross derivatives let us choose two coordinates and denote by 9i,52 the corresponding 
derivatives and by Jq^, Jq^ the corresponding integral operator as per the table (2-25): then we have 



d^d2:F = di / YdX 



(aiY)xdX 



(2-26) 



82 J di 



The important and conclusive remark now is that the order of the action of the integral operators appearing 
in the list (2-25) on Q. is immaterial because the kernel Q. is symmetric and -more importantly- because its 
residue on the diagonal is zero. This means that in exchanging two integral operators one may in fact acquire 
the integral of a total differential which is going to cancel either by integration or against the regularization. 
To illustrate the point we make two examples. 

Example 2.2 Consider two coordinates UK,a,vj^^: then 

d^K.c^d^j.p^ = Vj.fiUii.c^ ■ (2-27) 

In this case the two integral operators involve either residues (for K > 0) or (regularized) integrals. Either way the contours do 
not intersect and the double integral is independent of the order. 

Example 2.3 Consider the derivatives dug a ^""^ ^u^ai ^^^^ integral operators do involve intersecting contours, 

hence care must be exercised 

1 /■ „ .-K, 



2iKn 



■(C) 



f^(C,0 



The inner integral in fact does not need any regularization, so we have 

1 



(X-Qo 



'(C) 



f^(C,0 = 



= lim 



e->9c« 2iK-R 



(X-Qa 



1 



2iK-K 



-^(C)n(C,C)-;^(X(e)-Q. 



(2-28) 

(2-29) 
(2-30) 
(2-31) 



(The exchange of the order of the integrals gives a —2in5 supported at the intersection of the contours of integration). 
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Theorem 2.1 The free energy is given by the formula (we set uxfl '■= ux,oo, vjfl '■= vj^^d, f^o,o = ^o,( 
^0,0 = ^0,00 := for uniformity in the formute) 



rfi.a d2,a g 

= ^ UK,aUK,a + X] '"J:»^J>» + + X] ^J^J + ' 

a=OK=0 0=0 J=0 j=l 



- V resY^XdX 
2 ^ c 

- y resX^YdY 



(2-32) 



where 

Px := {oox,ga,0, a = 1,. . . ; j = 1, . . .} (2-33) 

Py := {ooY,Pa,r/i, a= 1,...; j = 1,...} (2-34) 

(see definitions of the properties of the points appearing here at the beginning of Sect. 2) 

Proof. First of all note that the expression is symmetric in the roles of X, Y after integration by parts and 
moving the residues to the other poles 

^ res Y^XdX = Yj res Y^XdX = ^ Yj resX^YdY , (2-35) 
Ce^x Ce»Y CgI'y 

where we have used that Px U Vy is the set of all poles of the differential Y^XdX. Now, the proposed 
expression is nothing but 

IT = ^^nK,a^^i^,a(YdX)-^ J^t;j,„Vj,a(XdY) + tT(YdX)+ (2-36) 

Q!=0_fS'=0 a=OJ=0 
9 1 

+ Y ( YdX) + - ^ res Y^XdX - t Y ^o,a (2-37) 

Suppose we compute a derivative w.r.t. u/j ^3: using the list of differentials (2-20, 2-21, 2-22) and moving the 
computation of residues over to Py -for convenience- before the differentiation, we have 



2du^,,T = UR,p(YdX) + YY''K,aUK,a(MR,f3^)+ (2-38) 

a=OK=0 

d2.c g 

YYl ^J,c.Vj,a {-UR^^n) + tT (UR^^n) + Y ^j^J {Ur,p^) - Y ^f^^nA^) = (2-39) 

a=OJ=0 i=i CG^y 

Y Y ^K,aUK,a i^) + YYl '^J^'^^M i^) + i^) + Y ^^^i " resXY(Q-|0) 

a=OK=0 a=OJ=0 j=l CS^^y ^ 



Note that the operator Ur^p involves residues at one of the points of Dx and hence commutes with the other 
residues when acting on the (singular) kernel also for the last term involving residues at Vy ■ 

^Thc set 2?x is the support of the pole-divisor of Y less the point ooy, and viceversa for 
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From the properties of il. and the definitions of the integral operators it follows that the differential acted 
upon by Ur^/^ is precisely YdX, namely 

YdX =Y.Y. ''K,MK,a {^) + Y.Y. + + £ ^i^i - . (2-41) 

a=QK=0 a=OJ=0 i=l C,eV^ ^ 

This can be seen by analyzing the singular behavior near the poles and the a-periods of both sides of the 
equality and verifying that they are the same*^ whence we have the desired conclusion of this part of the proof. 
The other derivatives are treated in completely parallel way. 

The derivatives w.r.t. Xj^Yj are a little different because there is no explicit dependence of JF from these 
coordinates. However this produces the correct result since, for example 

2dx,T = J^5^^^J,aVJ,,(-^^^^) + ^r(7^,^)+ (2-42) 

a=0 ^-=0 a=0 J=0 

9 

+ ^ejfj (7^^^) - J] resXY7e£(r?) = (2-43) 

= 7^^(YdX) , (2-44) 

which is consistent with our definitions (2-8). 

As a final case we compute the derivative w.r.t. Qa- here some care should be paid to the commutation 
of the derivative with the integral operators. Indeed dq^ does not commute with the integral operators 
Uk.oi^ K = 0, . . . but instead we have 

[Qq^MkA = = 1, . . . (2-45) 

{dq^M^A = . (2-46) 

While (2-45) is rather obvious from the definition of the integral operator, some explanation is necessary for 
(2-46). Expanding YdX in the local parameter at we have 

/■OOx 

YdX= -yi,,(X(e)) + co + ciz„ + 0(z2) + yi„(X(e)) (2-47) 



Therefore we have 



dQ^ YdX = dQ^ [\mi ^YdX + yi,,(X(e))^ = Sq^cq (2-48) 



Viceversa (recalling that dq^Za = —1) 

/■OOx /-OOx 



/■OOx /"OOx 

(9Q„Y)xdX = \im j (aQ„Y)xdX 



lim idq^co - ci + dq^cizc, + 0{zl)) = dq^co - ci (2-49) 

e — ^Qa 



^One should use that the behaviour near a pole of the last term in the LHS is, e.g. 

^ rcsXYf7~-d(Yy2',,(Y)) . 
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This shows that 



2m / X - (5o 



1,Q 



Using this and computing the derivative of T we obtain the desired result 

aQ„^ = resy{^(X)YdX . 



(2-50) 



(2-51) 



Finally, while the reasoning is mostly similar, the t derivative has an additional technical difficulty. First of al 
we have 



/•oox roox roox 

dt 4- Ydx = 4- 4- ^ + 1- 



J oc 



(2-52) 



OOy 



J OOy ooy 

The reason of the additional -1-1 is the fact that the local parameters near the two poles are different functions 
(here we set for brevity tx = t + UQ^a, tY = i + Y1 ^o,a) 



/■OOx 

dt-h YdX= lim dt 



YdX - (yi,oo(X) - txln(X))^ + {YVl^iY) - V2,oo{Y) - IyHY))^ 



lim 

e->ooY 
P — °°X 



P roox 



n + ln(X(/,)) - ln(Y(e)) + ( YV^(Y) - ^ ) (5tY)x 



€ J OOy 



= -j- -4- Q. + 1(2-53) 

J OOy OOy 

Moreover, whether we sum at the poles in Px or Vy, we need to interchange the order of the following 
residue/integral 



r>oox 

res XY / = lim 

e^oOY 



OOy 



OOy 



^2 oo("^)-"+E "0.a)Y"l + ... 
OOx ' II ' roox 

XYf] - Y(e)X(e) = -/- res XYJ] + t + V « 

^OOY „ ' 



res 

OOy 



Putting it all together we find 

2dtT = (T(YdX) +t-Y^ 7;o,a) + T(YdX) -t-^, vo,a = 2fi. 



(2-54) 



(2-55) 



The other derivatives w.r.t. to the moduli vj^a^Yj, Pa are computed in similar way by first rewriting the 
expression for T equivalently in the symmetric way w.r.t. the exchange of roles of X, Y. Q.E.D. 



Corollary 2.1 The Free energy satisfies the following scaling constraints 



l<a</3 



2T = Vx^ + ^ Uo^aUo,l3 + t ^ tio,a + 



(2-56) 
(2-57) 



l<Q</3 



where 



duK a 

a>OK>0 J=l 



dv 



J.oo 



a>l 



J=0 



dvja 



+ 
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a>OJ>0 ^^-^'^ K=l a>l \ K=0 OUK,a J 

Note that these formula give other representations of the free energy in terms of its first derivatives 
defined independently in (2-8). Moreover any convex linear combination will give another representation. 

Proof. The formulas can be obtained by explicitly computing the residues of Y^XdX at the various points or 
by the following straightforward argument. Consider the new functions X := X and Y = e^Y: the new free 

energy !F will be given by the same formula (2-32) in terms of the new objects. Taking ^ gives the first 

c=0 

formula. Some particular care has to be paid to the regularizations which involve subtraction of logarithms. 
The second formula is obtained in a symmetric way. Q.E.D. 

If we denote by Jq the integral operator associated to a derivative d, the formulas for the second order 
derivatives are written concisely 



9i92^= / / n + 59,,tSd,,d. (2-59) 

J di J 82 

In other words^ the Bergman kernel is the universal kernel for computing the second derivatives of the free 
energy and hence the two-point correlation functions of the matrix model in the planar limit. 

The third order correlation functions were computed in [2] for the case of polynomial potentials: since the 
reasoning is identical we only report the result. The key ingredient there is the formula that allows you to find 
the variation of the Bergman kernel under infinitesimal change of the deformation parameters. The formulas 
can be summarized as follows 

^ ^""^^'^^ " ip,aVc=t. dY(C)dX(C) - ^c=l dY(C)dX(C) 

™ ^ /■ V- o(e,c)^(p,c)^(r?,c) m.o^aiQMri.o 

id^Mi,,) = /Xffl. dY(C)dX(c) =^m. dY(c)dX(C) ^'-''^ 



k 



where Xk and yk denote -respectively- all the critical points of X and Y other than ooy,oox (namely 
dX(xfc) = 0, dY{yj) = 0). These formulas follow from Rauch variational formula [21, 15]. Note that dYdX 
in the denominator has simple poles at the hence the residues at these points do not contribute to 

the sum except for the cases where uq has a (double) pole at one of those points, namely only for the cases 

d = dx,, 

The final formulae for the third derivatives are simpler if we introduce the two kernels 



^^x(Ci,C2,C3) := -E/.- dY(C)dX(C) ^'-''^ 



k 



^Wc could dispose of the last term (enters only in d^J^) by subtracting ^i^; this would change the t-derivativc fi —> fj, -\- 1 making 
the formula for the first derivatives slightly different. Note that this does not affect the derivatives of order 3 and higher. 
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k 

This way one obtains 



^^v(Ci,C2,C3) .= >^^res 5Y(C)dX(0 ^'"''^ 



„ , 5^ = / Uk,JAjM$ > „ 5... . 9^ = / VA',a Vj,;3l^?) (2-63) 
Jd Jd 

du^ M^ = UK,aTT^^^ ■ dtdtT = Vj,aTTQ'}^ . (2-64) 
For all other third order derivatives one can use either kernels: 

(3) _ /■ f f n(3) 



SiSs^s^ = / / ^'4- / / ^x' (2-65) 

It should be clear to the reader that these formulas translate to residue formulas in the spirit of [19, 10]. For 
exannple 

d,^d,^d,,r = - V res ^^^^ = V res (2-66) 
^ C,=^k dYdX ^ c=yk dYdX 

k k 

We should stress, however, that although the present moduli space can be embedded as a submanifold of the 
moduli space considered in [19], the coordinates {uk,vj) that are relevant to the matrix-model applications 
are of a different sort and -resultingly- the free energy it is not the same function as the tau function of the 
Whitham hierarchy in [19]. 

3 Residue formulas for higher derivatives: extended Rauch- variational 
formulae 

It is clear from the previous review of the material that in order to compute any further variation we must be 
able to find the variation of the kernels ily^ and this step will produce three kernels 

r^YY > ^YX ' ^XX ' (3"-'-) 

according to which variable Y or X we keep fixed under the new variation. The reason of this plethora is 
essentially that the variations of the basic differentials are performed more easily either at Y or X fixed: 
for instance if we compute the variation of lok = Uk{^) at X-fixed we obtain 

{duK^MO = res — (aO)x = -\^ res (3-2) 



K ' ^c=a.. dY(C)dX(C) 



k 

whereas 

{duJK,oo)Y = res X^-1(5X)yJ^ + -^{dn)Y (3-3) 

This is in fact a manifestation of the thermodynamic identity for differentials 
Lemma 3.1 The variation of a differential at Y and X fixed are related by tiie following formula 

(dio)^ = (aa.)x + d (^(aX)^) = (9c.)x - d {^y) (3-4) 
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Proof. Writing uj = /dY = (/dX we have 



Y 



(55)x + ^(5X)y 



dX + 5d(aX)Y = {du)^ + d (5(5X)y) . (3-5) 



Since g = w/dX and (9X)y = — i^a/dY we have the assertion. Q.E.D. 

Using Lemma 3.1 and trading the residues at the Xk over to the others (at the y£,CiOOx) one can check 
directly that the formulae (3-2, 3-3) are consistent. 

It should also be clear that the variation of the numerators of x ^''^ obtained by simply applying the 
product rule and the previously listed appropriate Rauch formulas. The only new ingredient is the variation of 
the denominator of x explained below. 

Suppose we want to perform a variation 9 at X fixed of one of the two kernels; when we need to compute 
the variation of the denominator we need a formula for d-^. We should think of the expression as a 
meromorphic vector field on the Riemann surface and the variation is the vector field 

/ 1 \ d((OY)x) 1 ,(^d\ .... 

HdY = dY^ = "dY^^ (dX j 

Now, the differential of the function ^ can be expressed as a residue using -once more- the Bergman kernel 
according to the following 

Lemma 3.2 Let F be a (local) meromorphic function: then the differential dF can be obtained by 

dF(C) = resn(C,C)nC) . (3-7) 

The proof is very simple using a local parameter near the point ^ and the asymptotic expansion of the Bergman 
kernel. 

Combining Lemma 3.2 with (3-6) we have the new variational formula 



dY)x 

1 

dX 



dY2(e)c=€ ^^'"MX(C) 
1 resO(C,0^S^ (3-8) 



dx2(0c=« '^'^My(C) 



where the different sign in the second formula is due to the fact that (5X)y = — cjg/dY. 
Let us summarize the rules of the calculus: 

1. The variations of any differential can be performed at X or Y fixed, the two being related by Lemma 
3.1. 

2. The variations at X-fixed of the vector fields 1/dY and viceversa are given by eq. (3-8). 

3. The variations of the Bergman bidifferential Q are given by eqs. (2-60) 

The choice of variable to be kept fixed Y vs. X is ultimately immaterial. However formulas can take on a 
significantly more involved form if one chooses the "wrong" way of differentiation. We are going to practice 
this calculus and compute the fourth order derivatives explicitly. This will also provide us with relevant formulae 
for the four point correlators of the planar limit of the two-matrix model. 
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3.1 Fourth order 



To illustrate the method we compute the fourth derivatives w.r.t. UK,a, '>J'L,i3, um,"/, un,5 (which we will denote 
in short hand by subscripts MiN iL iK only). We start from the expression for the third derivative 



5=x. dYdX 



(3-9) 



It is quite obvious from the considerations around Eq. (3-2) that the extra derivative is most easily computed 
at X-fixed: 



dK^ALN^L 



res 



i=x, dYdX 



dY 



Using now Lemma 3.2 and the variational formulae (2-60) we obtain 



M,N,L 



The computation could end here, since we have successfully expressed the derivatives in terms of residues of 
known differentials: however this expression is not obviously symmetric in the exchange of the indices, whereas 
it should be since it expresses the fourth derivatives of the free energy. The expression is symmetric, but not 
at first sight. In the double sum the order of the residues is immaterial only for the non-diagonal part: for the 
diagonal part of the sum the residue w.r.t. ( must be evaluated first. The non-diagonal part of the sum is 



'^L(6^Af(0 o^e A\^Af(C)'^i^(C) , I fA,r at\ 



^€=x.c=-. dY(OdX(0 



'dY(C)dX(C) 



(3-13) 



where the order of the residues is -as we said- immaterial because they are taken at different points. This 
term corresponds diagrammatically to 




® ® 
and is manifestly symmetric in K, L, M, N. The diagonal part is not manifestly symmetric but in fact we 

are going to show that it is. The diagonal part of the sum together with last term is made of the following 

residues 



A- dY(Odx(0 dY(c)dx(c) ■ + - ^ - ^) + "^dY(eFdx"(c;'' m ""dx(0 ^^-^'^ 

where we stress that the residues w.r.t. ( have to be evaluated first. For instance, a rather long computation 
in the local coordinate z = -y/X — X(xa,.) gives 

ILMN K" + L M N" K + L M" N K + L" M N K + L M N K Sb IKLMN Y'" 



(Y') 



l\2 



(Y') 



/^3 



(3-15) 
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where the short-hand notation is as follows 

lul = L(z)dz , ujk = K{z)dz , ujm = M{z)dz , ujn = N{z)dz (3-16) 

and 5^(2;, 2:) is the projective connection, and all quantities are evaluated at 2; = 0. 

3.1.1 4-point correlator 
This is the formal expression for 

5^ J' 

0^1(^1)^1(92)^1(93)^1(^4) 

where the formal operator S/6Vi{q) is defined by 

^ = V q-^-^K-^ . (3-19) 

5Vi(g) duK,oo ^ 

By summing the four indices of the above derivatives (at least formally) we obtain 

Rffl{qi,q2,q3,q4.)dqidq2dq3dqi = J^xx(C(9i), ((^2), ((^s), ((^4)) (3-20) 
where ({q) is the solution of X(C) = g on the physical sheet of the cover X : ^ CP^ and 

J7(^)(1 2 3 4) - Trc. ^(^'g)"(^'^)^(3,e) ^.^^^. .0(C,4) 

iZxx(i,A3,4j - 2^res^ dY2(e)dX(e) C=l dX(C) + ^ ^ 

+ V res V res ^^MlilMlf^f. ^(^- ^ + (1 ^ 2) + (1 ^ 3) (3 22) 
4-/=t ^C=t dY(C)dX(C) ^^^'^'^^ dY(OdX(e) + (1 ^ 2) + (1 ^ 3) (3-22) 

Note that this kernel is symmetric in the four variables although not at first sight, but by the same consider- 
ations as before. In a similar way one can obtain the other four point correlator 

-RK(pi,P2,P3,P4) := ... . u//u!^^ ^TA^ ^ ' (^"^3) 

0V2{P1)V2[P2)V2{P3)V2{P4) 



where S/6V2{p) is defined similarly as before by 

The derivation of the formula is completely parallel hence we only give the final result 

^0^1(^1 'P2,P3,P4) dpi dp2dp3dp4 = f^YY(C(Pl),C(P2),C(P3),C(P4)) (3-25) 

where ^(p) is the solution of Y(^) = p on the physical sheet of the cover Y : — > CP^ and 

o(^)(i 2 3 4) - y-,,, »(i^O^(2,OT3,0 ^.^^^. ..^(C,4) 

iVY(i>A3,4j - ^res ^Y(e)dX2(e) C=| ^^'^^dY(C) + ^ ' 

^ ^ ^^(1, 0^(4,0 ^,, ^, 1^(2,0^^(3,0 ^ oWon^N 

+ E/- E^Lt dY(OdX(0 "^^'^^ dY(OdX(0 + (1 - 2) + (1 ^ 3) (3-27) 
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3.2 "Mixed" fourth order derivatives 



As a further example we compute the derivatives w.r.t. ul,um,vn,vk- we leave the derivative w.r.t. vk last 
and perform it at X-fixed 



res 



UJLUJMUJjy 



E 



res 



res 



dYdX 



dYdX 



dY 



(3- 



dYdX dYdX VdYdX; dY2dX VdXX^ 



E 



res 



wa./(C)'^l(C) 



C=x. dY(C)dX(C) ^?=?/. 



E 



res 



E 



res 



dY(OdX(0 



+ (L ^ M) + 



+ > res 



C=x. VdY(C)dX(C) ^^=x, dY(C)dX(C) 

dY(C)dX(C) VdY(C)dX(C)y' dY(C)2dX(C) VdX(C)y 



(3-31) 
(3-32) 
(3-33) 



Note that in order to compute effectively the derivative of tOj^ at X-fixed we have used Lemma 3.1. Once 
more one can check that the resulting expression is symmetric \n N ^ K and M ^ L. The only terms which 
do not have this symmetry at first sight are the diagonal part of the double sum over the Xk together with 
the last term: 



LMNK 



Eres > 



res 



WA/(C)t^L(C) o^,^ ^^^iv(0'^A■(0 



-E 

-E 



S=y, dY(C)dX(C) 



■ppc -ppc 

C=x.«=x. VdY(C)dX(C) 



dY(C)dX(C) 



res 



dY(e)dx(e) 

wa/(C)wl(C) , f ^M'^kiO 



(=x, dY(C)dX(C) VdY(C)dX(C) 



+ 



res res 



C=x. V?=^. dY(C)dX(C) 



dY(OdX(0 dY(C)2dX(C) VdX(C) 



(3-34) 
(3-35) 

(3-36) 
3}37) 



The same considerations about symmetry done previously apply to the sum on line (3-37) as well. 



3.2.1 4-point correlator 

Using the above computation we can compute the following four-point correlator 



-^2,2(^1' 92, Pi, P2) : = 



Performing the multiple summation we find 



SVi{qi)Vi{q2)V2{pi)V2ip2) 



-R2^2(9i,?2,Pi,P2)dgidg2dpidp2 = ^^xy(C(9i), C(92), ^(pi), ^(^2)) 



(3-38) 



(3-39) 
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where ^{p) is the solution on the physical sheet of Y(^) =p and 

fi(^) (1 2 T 2) - V res V res ^^MMMlf^^. O^lMMlil + (3-40) 
U^^(l,2,l,2) - 2^^«;s^2^ res dY(C)dX(C) ^^'^^ dY(OdX(C) + 

^C=2:.e=a:. dY(C)dX(C) dY(OdX(0 

+ E dY(Odx(0 dY(Odx(0 + ('-''^ 
- L I^.Lt dY(Odx(0 ^) dY(Odx(0 + dY^(Odx(0 I, dx(0 7f ^ 

Repeating the derivation from the beginning one can realize that there is no need of any other kernel for 

6^ J' 

R^} (91 , Q'2 , gs , P2 ) : = .T. / NT./ f f \ ' (^"'^'^) 
oVi{qi)Vi{q2)Vi[p3)V2[pi) 

which is given by 

^3^1(91, 92, 93, Pi)dgidg2dg3dpi = l^xx(C(9i),C(92),C(93),C(pi)) ■ (3-45) 
3.2.2 Summary of all fourth derivatives 

These three kernels are sufficient for us to write all fourth derivatives compactly as some new residue formulae 
(note: the order in which the integral operators appear is to mean that they are applied to the variable that 
appear in the corresponding position in the kernel) 

du^ du, d,, d,,, T = UkUjVlVm^'^y (3-46) 
du^ d.,, d^^T = UKTTV.M^-y (3-47) 

du^did2dsJ' = UK [ [ [ 1^5x (3-48) 



d,,Ad2d3T = Vj [ [ [ (3-49) 

J 81 J 82 J 83 

did2d^d,T = [ [ [ I n4^-y= [ [ I / i^S^ (3-50) 



'81 J 82 J 83 J 8i J 8\ J 82 J 83 J 8i 

where the symbols here mean derivatives with respect to variables not included in the previous items of the 
list. 

3.3 Higher order correlators 

The computation of any derivative of any order is just a matter of application of the "rules of calculus" 
outlined previously; in this fashion one could obtain residue formulae for any derivative and possibly develop 
some diagrammatic rules to help in the computation. We leave this exercise to the reader who may need it 
for his/her application to a specific problem. The formal "puncture" operators 
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act as follows on each term 

'^^^^^<5yi(X(l)) VdY(2) j " dY^'^' V dX(2y J " dY^ |=l dX(0 ^ ' 

s f ^ \ f m,2) \ i_ n{2,on{^,i) 

'^^^^'6V2{Y{l))\dX{2)) ~ dY2(2)'^'UY(2)y dX2(2) P dY(e) ^ ' 

Combining these "rules" it is easy to obtain any correlator: the resulting expression will be symmetric in the 
exchange of the variables, although to recognize this some careful analysis of the residues is required. 

3.4 The Equilibrium correlators 

The derivation of the multiple derivatives of the equilibrium free energy Q follows the same lines and the 
results are the same formulae with Q replaced by Q (clearly there are no derivatives w.r.t. the filling fractions 
ej which are now dependent functions). In general the rules of calculus for Q are the same as the rule of 
calculus for T with all the instances of the Bergman kernel replaced by the dual kernel il. 

A Explicit form of the regularized integrals 

In this section we provide explicit formulae for the regularized integrals used in the definition of the Free energy 
and the r function of the previous section. 

The main tools are the following properties which were used in the proof of the derivatives of the free 
energy 

YdX = ^K,aUK,c. i.^) + Y.Y. ^J^^yj^o. {^) + tT (0) + ^ e^e^ (SI) - res XYO (A-1) 

a=ox=o Q=oj=o i=i Ce^Y 

XdY = -Y^ll ^K,aUK,a {^) - ^J,c.Vj,a i^) " tT (0) - ^ e^Sj (fl) - res XY^-2) 

a=OK=0 a=OJ=0 j=l (eV^ 

Let us compute ^^'^ YdX according to the original definition of regularization: since the operator -J^^ 
commutes with the integral operators/regularizations in (A-1) we obtain immediately 



4- YdX = n,,,«Wi^,„ / +5^J^fj,aVj,„ / n\ + 

•'1" a=OK=0 ^-'1" ^ a=OJ=0 ^"^^o / 

(/■oox\ 9 / rooxX /-OOx 

/ Q \ + "Yej£j i ^) - Yl resXY / Q (A-3) 

The differential is the unique normalized differential of the third kind with simple poles at (?a,oox and 
residues -respectively — hi, —1. To simplify formulae let us define for any two points ^, i] the following function 



C fV 

n 



Co 



a = (A-4) 
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This is a multivalued function around the 6-cycles; on the simply connected domain obtained by dissection of 
our surface, A^^^ has a simple pole at ^ and a simple zero at r/. It is defined up to a multiplicative constant 
(depending on the base-point for the outer integration), which however will not affect our result. With this 
definition we have {qq := oox) Po ■= ooy) 

duo^:F=4- YdX = - VresFi,5(X)— + Vres(y2,/3(Y)-XY)— + 



where we have set A := Ag^ oox ^i^d 



ln7oox := In(Ag^ ooxX) 

ln7,„ := lim ln(A,„,oox(X-Qa)) (A-6) 

The formulas for the derivatives w.r.t. vq^o. are obtained by interchanging all the roles of X, cxoxjQa with 

Y, OOY, Pa- 

Finally the formula for the t-derivative 

/■OCX A\ A\ 

dtT = 4- YdX- V^;o.a = - Vresyi,5(X)— + Vres(y2,^(Y)-XY)— + 
Joo^ „ A f-' A 



a 5=0 /3=0 



+ 

where -this time- 



A . AqoyjOOx 

ln(7oox) := lim In (AX) 

e^oox 

ln(7ooY) := lim In ( - ) (A- 



The extra " ^Q,uo,a " which cancels with the same term in the expression for /i is due to a careful analysis of 
the regularization prescription for the following term in the computation 



res 

OOy 



/■OOx / /"OOx \ /-OOx 

XY / Jl= lim / resXYO + X(e)Y(e) =-t-y7;oa + -/- res XYJ] . (A-9) 
Note that, in all these formulae, the 6-periods of ^ are the Abel map of the two poles of this differential. 

B Example: one cut case (genus 0) and conformal maps 

The formulae for the derivatives simplify drastically in case the curve is a rational curve. In this case, 
introducing a global coordinate A (as explained in [1, 2]) with a zero at ooy and a pole at oox and suitably 
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normalized one can always write the two functions X, Y as 

X = 7A + I; ^i,,ooA-^ + E E - Aa)-^-^ + E 

K=0 a K=Q i=l -^'^ 

Y = T + E ^i.-A^- + E E - ^-y'~^ + E Y^:^ (B-i) 

J=0 a J=0 i=l •''^ 

Q„ := X(A,) , Ps := Y(A5) . (B-2) 

The parameters ^,Aj,Fj and Bj,Gj,j = l...,s are not independent but are constrained by the following 
set of linear equations (in ^, Aj, Bj, Fj,Gj) 

X'(A,,x) = 0, j = l,...,r ; Y'(A,- y) = , j = 1, . . . , s . (B-3) 

As we have already mentioned, the equivalent of the Bergman kernel is simply 

dXdfj, 
(A^2 

This is the kernel of the derivative followed by projection to the principal part. For example the differentials 



^(A,/^) = 7f^- (B-4) 



Wif,oo(A) = - res ^^^J7(A,m) = ^(X^), dA, (B-5) 
^=00 A A 

where the it-subscripts mean the polynomial or Laurent part of the expression enclosed in the brackets. Similar 
completely explicit formulas for all other differentials in the lists (2-20, 2-21, 2-22) are left to the reader. 

The coordinates are given by the usual formulas (2-3). The free energy can be written in a quite explicit 
form using the following simplifications due to the existence of a global coordinate A (Aq := 00, Aq := 0) 

duo,.:F = -E/_^s ^i(X)y^ + E/_T(^2(Y)-XY) 



A=A^ A — Aq _ A=A^ A — Ac 

/3>0 /3>0 



+ E ^o./3^"(7(A/3 - Aq)) -t^o^aln f 7^ ) ) + E ^O'/^ ^'^ ( \ , \ + tln{Xaj) (B-6) 

/3^{5,0} 



(A^-A5)7y ' V 7 ^ ^ ' VA/3-A5 

A 



since Ag^ = x-x ^"^"^ Ap^^ooy = x-x- Moreover, using this time Aoox.ooy = A and formula (A-7) 



+ ^ no,/3 In (A;37) " E ^0-/3 ( " ) + * (^') + * (B"^) 
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By computing the other residues one can get explicit formulas for the Free energy in terms of the uniformization 
(B-1) and using Thm. 2.1. 

Denoting as before by Xk and y£ the critical points of the functions X and Y^*^ respectively we have as 
example of fourth point correlators 



+E 



(Ail - Xk) (/i2 - Xk)' 



1 



(Ai3 



Y>{xk)X."{xk) {Xk - Xrf 



6Y'"(X")^ 



2Y'(X'")^ + 2Y"X"X"' 



Y'(rE,)X"(x.) 
- 3(X")^Y"' - 3X(^^)X"Y') 



(/il - X/c)2(^2 - 2;fc)2(/i3 - Xfe)2(^4 - Xfc)^ 



+ (/ii ^ /is) + (/ii ^ //4) +(B-9) 

(B-10) 



+ 



18(X")2Y'((/.i-Xfc)-2 + cyc) 



2X'"X"Y'((/ii-Xfc)-i + cyc) 



(/il - Xfc)2(^2 - 2;fc)2(/i3 - Xfc)2(^4 - Xfc)^ (/ii - Xfc)2(^2 " 2;fc)^(^3 " 2;fc)2(//4 - Xfc)^ 



(B-11) 



Here the expression looks more complicated than necessary because the derivatives are taken w.r.t. A. 

Higher order correlators are of increasingly cumbersome expression but in principle they are easily computed 
using the general calculus outlined in the main text. 

B.l Conformal maps 

A further simplification of the formulae arises in the case the functions Y and X above describe the Riemann 
uniformization and its Schwartz reflected of a simply-connected domain V in the X plane []. We recall that all 
our formulas can be easily adapted to the description of simply and multiply connected domains (the number 
of connected components being the genus of the curve) by taking the curve as an M-curve in the sense 
of Harnack [14]: namely a curve with an anti-holomorphic involution ip : T,g ^ T,g having g+1 contours of 
fixed points and such that 

X(C) = YMC)). (B-12) 

In genus zero and with the normalization used in the previous paragraph for the uniformizing coordinate, the 
anti-holomorphic involution would be A — > =. The two functions Y and X then satisfy 



X(A)=Y( = 



(B-13) 



Since X(A) is now the uniformizing map of a simply connected domain V it follows from the general properties 
of such maps that X maps biholomorphically the outer region C\V to the outside of the unit disk in the 
A-plane. This means that the zeroes of dX all lie inside the unit disk and hence the zeroes of dY (which is 
the Schwartz function of the domain) all lie outside. 

The Free energy of the two-matrix model under this reduction vk = uk, reduces to the tau-function of Jordan 
curves studied by Zabrodin at al. [17, 18, 20, 23, 24, 25] as explained in [1, 2] 



1 

4^ 




d^Xd^Xln 



V JV 



1 

X 



X 



(B-14) 



'^"Note that the set of points {Aj_x} is a subset of the {xj;} and similarly for the {Aj y} and {yk}- 
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The coordinates uj^ are the so-called exterior harmonic moments of the region 




— / X^^'dX A dX . 
2z7r J J) 



(B-15) 



(B-16) 



The Free energy is in this case a real analytic function of the harmonic moments T = {uK,UK,t) and the 
previous formulas for the fourth derivatives^^ can be translated into contour integral-formulae which in turn 
could be written in terms of the Green's function of the Laplacian for the given region. It is also clear that 
effective formulae can be obtained for the multiply connected domains which correspond to higher genus 
-curves considered in this context. 

C An extended moduli space 

The moduli space considered in this paper could be easily extended in the spirit of [19] by considering instead 
of functions X, Y some normalized second-kind differentials dX,dY: this generalization has probably no 
relevance in the context of matrix models, nevertheless we sketch the main extra features. The practical 
difference is that now we may still think of multivalued functions X, Y with the properties 



whereas the functions have no multivaluedness along the a-cycles. The rest of the description of the moduli 
space is exactly as in Sect. 2. Note that this moduli space is "larger" than the moduli space of [19] because 
we are also considering the position of some zeroes of our primary differentials. 

After dissection of the surface Hg along the chosen cycles {aj,6j}j=i^ ^ and along the fixed contours 
between the non hard-edge poles we obtain a simply connected domain over which we will consider the 
functions X = J dX, Y = / dY. In this domain the same asymptotics as in (2-3) are valid (where the 
"potentials" are discontinuous across the cuts along which we have dissected the surface). The free energy 
(we should probably call it rather "tau" function) would be defined by the same formulas (2-8) except for the 
fact that the e^-derivatives should be replaced by the formulas below and we should consider the derivatives 
w.r.t. the extra moduli Aj,Bj 



Y(C + bj) = Y(C) + Bj , 



(C-1) 
(C-2) 




(C-3) 



(C-4) 



(C-5) 



^Thc third derivatives were computed in [25]. 
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The equivalence of the two last lines is given by integration by parts. Also, the last integrals may seem to 
depend on the base-point of integration: in fact they do not as one may check by computing the differential 
ate. 

Besides the differentials considered in (2-20, 2-21, 2-22) one also has 

{dB,Y)^dX = ^ ([ Xn =: B,in) (C-6) 

These formulas are obtained by noticing that (5B^Y)xdX is a holomorphic multivalued differential with 
monodromy only around the corresponding 6-cycle 



(5B,Y)xdX 



= -<5,fcdX(C) . (C-7) 

c 



The integral formula has the same properties and hence we have the equality. The reasoning for (3^^X)YdY 
is symmetric. Note that -using the thermodynamic identity- we have 

{Ob, X) YdY = - ^ ^ XO . (C-8) 

The considerations to prove the compatibility of the above equations are similar to the previous case with 
one notable exception we want to bring to the attention of the reader; in the computations of the second 
derivatives one is lead to considering integrals of the form 

XS Yn , J) xJ) Xn , i Y i YQ. . (C-9) 
These integrals do not depend on the order only if j ^ k: in fact we have 



X(b Yn = (t Y (b Xn + 2iTT d) XdY 

xixn = (b X ([ xn + 2iTT ([ xdx = ([ xJ) xn 

I Y l YVL = I Y l YVL + 2m I YdY = / Y / YQ (C-10) 

J aj J aj J aj J aj J aj J aj J aj 

Another kind of integrals that one encounters are of the type 



Y * n = 2iTT 6 YiUk (C-11) 

aj J bk J aj 

Here one has to use the following rule for exchanging the order of the integrals: suppose that a specific choice 
of the homology representatives of aj and hj intersect at the point ("o. then 

/ F(c)/ nic,o=(f (F(c)-F(Co))/ n{c,0 + F{Co)<f i ^((,0= (C-12) 

= 1 I {F{C)-F{Co)MC,0-2i7r6,kF{Co) = 2i7rSjkF{Co)+ (f I F{C)n{C,^) (C-13) 
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By following similar arguments one can prove that 

9 

2T = 2To + ^(Aj-Aj + BjMj) (C-14) 
i=i 

where To is given by the same formula (2-32) (with the new meaning of Tj, though). The proof rests on the 
identity 

YdX = 5^5^ni^,<,Z^x,a(f^) + 5Z5Zi^j,aVj,„(0)+ (C-15) 

0=0 K=0 a=0 J=0 

+tT{n) + y lejSj {Q) + — J) X.Q + — J) Yn] + V resXYO , (C-16) 

which is proved as before by matching the singular behaviors of both sides at all possible singularities and by 
checking that both sides have the same multivaluedness around the a and 6-cycles and the same periods. 

C.l Higher order derivatives 

In order to write compactly the second derivatives let us denote by d any derivative w.r.t. one of the parameters 

UK,a,Qa, Xj,'^j,aT Pai^j- Beside the second derivatives already computed, the new ones are given by the 
formulae 

s 

dA,dB,T = AjBkn + ^ek 

dAA.^ = -^j^k^ + 

dB,dF = Bj [ n ■ dA,dT = Aj [ n (C-17) 

' Jd ' J d 

We remark that the order of the integral operators acting on Q. is relevant because $7 is singular on the 
diagonal: for instance 

A^Bk^l = BkA.Vl - . (C-18) 

In order to compute all higher derivatives and loop correlators we need to specify the relevant additional Rauch 
variational formulae: besides those considered in (2-60) we need the ones related to the extra moduli 



(9b,J^)x(1,2) = ^jf^ X(e)J7g)(l,2,e) 



dY(0 

(9^^.0)y(1,2) = Y(e)0§^(l,2,0 (C-19) 
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We briefly justify these formulas. Suppose uj is any of our differentials and consider the function oj/dX (or 
symmetric argument for Y). This function has poles at the zeroes of dX and possibly constant monodromy 
around a 6-cycle. Thinking of it as a function of X the monodromy condition reads (c is or 1 depending on 
the case chosen, but the argument is unaffected) 

^(X + yl,)-^(X) = c (C-20) 



Taking the derivative w.r.t. Aj at X-fixed we have 

(<9A,t^)x 



c 

Considering with some care the singularities at the zeroes of dX and this multivaluedness one gets 

id.uMn- V res ^^^^Hii^lMM + ^ / ^{C.^MO ,^ 
(dA^^m - - ^ res^ dX(OdY(0 + ^ t, dX(0 " (^"''^ 

This gives the previous extended Rauch formulas. 

Using these expressions for the variation of the Bergman kernel one can obtain all third derivatives. Besides 
those already considered in (2-63, 2-64, 2-65) we find also 

dA, dA, dA,T = AjAkAen^ ; Ob, Ob, Ob.T = B^BuB,^'^ 

(3) 
X 



dA,dA,dF= [ AjAk^^^^ ; dB,dB,dT= [ BjBk^ 

Jd Jd 
n T- A t2 c r.{3) , 1 / Bk{^)£j{^) Sjk6kl 

dAMd.r = + — £^ dX - 

2^7^/„^. dY 

F> T Ann o(3) ^ 1 / Bu^UkA^) 
dA^dB^du^^^T = AjBkUK,a^x + 2i7r J 



dX 

dA,du^,^duj,,T = AjUK,aKj,i3^'x. + ^ f dX 

dB^dy,_A'.,,T = BjVK,aVj,pn'Y + ^ dY 

dAAd^= I TAjQ^^^ ; dB,dtd:F= [ TB^^^^ . (C-23) 

Jd Jd 

C.1.1 Order four and higher 

it is clear that the formulae become rather long due to many case-distinctions. However the reader should be 
able to compute any derivative of order four or higher by using the same rules of calculus outlined in the main 
text, with the additional Rauch formulae (C-19). 
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D General definition of regularized integrals 



Let a; be a meromorphic differential with poles at the points (p, p = 0, . . .. Let Zp be chosen and fixed local 
parameters at (p. Let loj be the Abelian differentials of the first kind nornnalized w.r.t. the a-cycles of a given 
choice of basis {aj,bj} in the homology of the curve. Then we have 

LO = 

p>OK>l^' p>l ^ -' j=l 

The regularized integral from to is defined for a homology class of contours in the punctured surface: 
in general one has to dissect the surface along the a, ^-cycles and along a set of mutually non-intersecting 
segments joining the poles of m in such a was as to have a simply connected domain. Choosing an arbitrary 
path within this simply connected region and joining the two chosen points we have (supposing that both ^,ry 
are poles of uj) 



n>^)K>^ n>^ ^ n=^ V"-^ J 



p>ois:>i 



A 




(D-2) 



7? := limln{^) (D-3) 
7^ := limln(A(e)z^(e)) . (D-4) 

Some remarks are in order: the function ln(A) is defined as any antiderivative of the normalized third kind 
differential 17, which has residue —1 at ^ and residue +1 at r]. Hence A has a simple zero at ^ and a simple 
pole at rj (in the simply connected domain). Also, A is defined up to a multiplicative constant depending on 
the base-point of integration: the final formula for the regularized integral does not depend on this constant. 
A can be written explicitly in terms of a Theta function and the 6-periods of ^ are the difference of the 
Abel-map between ^ and r]. 

In the more general situation of the extended moduli space studied in Sect. C we had also some multival- 
uedness of the type 

u;{C + b,)-uj{0=dH,{0 (D-5) 

where dHj{Q, j = l,...g are meromorphic differential of the second kind with vanishing a-cycles. The 
formula for a regularized integral is easily adapted: the main observation is that (D-1) now needs on the 
R.H.S. the following extra term 

^ = (d-i) + ^E/ 

i=i -^"^ 

and consequently the formula for the regularized integral is 

1 f „ dA 



? j=i 



UJ 
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